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Abstract 

We  derive  and  implement  analytic  solutions  for  the  description  of  batteries  subject  to  cyclic  voltage  and  current  sources  of  varying  amplitudes, 
consistent  with  vehicle  drive  events.  An  equivalent  circuit  comprising  a  resistor  in  series  with  a  parallel  resistor-capacitor  combination  is  employed. 
The  method  we  develop  is  applied  to  a  high-power  lithium  ion  cell,  thereby  allowing  us  to  obtain  parameters  and  overall  characteristics  useful  for 
(1)  representing  the  battery  in  system  (e.g.,  vehicle)  models  and  (2)  providing  a  quantitative  means  of  comparing  and  classifying  battery  systems. 
Closed-form  analytic  solutions  derived  using  Fourier  series  are  compared  with  experimental  data. 

©  2007  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

To  facilitate  the  integration  of  traction  batteries  within  vehi¬ 
cle  systems,  procedures  must  be  developed  to  characterize  the 
devices  and  extract  relevant  parameters.  The  control  algorithms 
employed  in  hybrid  electric  vehicles  (HEVs)  rely  on  simplified 
models  that  are  linear  and  admit  rapid  calculation,  commen¬ 
surate  with  robust  vehicle  operation  and  embedded  controller 
capabilities  (execution  times  and  storage).  Such  control  algo¬ 
rithms  are  used  to  determine  the  state  of  charge,  state  of  health, 
and  power  capability  of  the  energy  storage  system  [1-5].  The 
United  States  Advanced  Battery  Consortium,  in  conjunction 
with  the  FreedomCAR  initiative,  has  provided  target  specifi¬ 
cations  for  HEV  batteries  and  associated  test  procedures  [6]. 
The  test  procedures,  however,  do  not  lend  themselves  to  the 
straightforward  determination  of  the  parameters  associated  with 
the  equivalent  circuit  depicted  in  Fig.  1,  which  is  often  used  to 
correlate  battery  performance  and  is  used  by  many  automakers 
(including  GM)  to  represent  battery  systems  in  full  vehicle  mod¬ 
eling.  The  parameters  R ,  Rc t,  and  C  are  expected  to  vary  with 
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temperature  and  SOC,  and  a  method  to  assess  such  variations  is 
needed. 

By  controlling  the  cell  potential  (and  enforcing  current  lim¬ 
its),  one  can  ensure  that  the  cell  state  is  well  controlled  during 
the  test  conditions.  Cyclic  voltammetry  [7],  wherein  the  cell 
potential  is  changed  linearly  with  time,  renders  an  easily  imple¬ 
mented  and  interpreted  procedure  for  analysis.  There  are  also 
times  when  electrochemical  systems  are  best  analyzed  using  a 
triangular  current  sweep  excitation  source  [7],  particularly  when 
one  wants  to  avoid  high  current  passage.  In  this  work  we  pro¬ 
vide  analytic  expressions  for  cyclic  voltammetry  and  triangular 
current  sweep  chronopotentiometry  of  electrochemical  cells  that 
can  be  represented  by  the  equivalent  circuit  depicted  in  Fig.  1. 
An  analogous  approach  was  recently  presented  by  the  authors 
for  supercapacitor  cells  [8].  In  that  case,  a  series  RC  circuit  was 
employed;  i.e.,  Rct  infinite  in  Fig.  1.  All  formulae  derived  in  this 
work  for  cyclic  voltammetry  transform  to  those  of  Ref.  [8]  when 
Rct  is  taken  to  be  infinite. 

When  developing  an  approximate  model  for  a  complicated 
system,  one  must  ask  “can  the  data  be  well  represented  by 
the  model  if  the  model  parameters  are  held  constant?”  If 
the  model  cannot  correlate  a  broad  range  of  data  with  suffi¬ 
cient  accuracy,  then  additional  complexity  must  be  incorporated 
within  the  model  representation.  We  include  a  number  of 
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Nomenclature 

C  capacitance  (F) 

F  Faraday’s  constant  (96,487  C  mol-1) 

I  current  (A) 

Q  charge  on  the  capacitor  (C) 

r  R/Rct 

R  high-frequency  resistance  (Q) 

Rct  Effective  interfacial  resistance  (Q) 

Rg  gas  constant  (8.314  J mol-1  K) 

SOC  percent  state  of  charge  (energy  content  in  the 
battery  relative  to  the  energy  content  upon  full 
charge) 
t  time  (s) 

T  temperature  (K) 

U°  defined  in  Eq.  (1) 

V  cell  voltage  (V) 

Vo  open-circuit  voltage  (V) 

Greek  letters 
a  RC/0 

v  frequency  of  excitation  source  (s-1) 

0  half-cycle  duration  (s) 

r  t!0 


experiment-theory  comparisons  over  a  broad  range  of  operating 
conditions  using  the  same  set  of  parameter  values  in  this  work. 
An  important  implication  associated  with  the  successful  com¬ 
parison  of  the  analytic  solutions  with  experimental  data  is  that 
control  algorithms  based  on  the  equivalent  circuit  of  Fig.  1  can 
be  expected  to  work  well  in  terms  of  correlating  (adaptively)  the 
performance  of  the  battery  system.  While  the  algorithms  can  be 
implemented  to  regress  adaptively  the  parameter  values  during 
vehicle  operation,  it  is  necessary  to  place  bounds  on  the  adapted 
parameters  for  robust  operation  [9-18],  and  the  base  parame¬ 
ter  values  are  also  required  to  assess  the  state  of  health  of  the 
battery.  In  a  related  vein,  we  note  that  several  researchers  have 
attempted  to  infer  the  SOC  by  first  correlating  the  equivalent 
circuit  parameter  values  as  a  function  of  SOC  and  then  measur¬ 
ing  the  parameter  values  (e.g.,  by  ac  impedance  methods)  of  the 
battery  in  different  states  [19-25]. 

The  intention  of  the  simplified  model  presented  in  this  work 
is  to  provide  rapid  parameter  identification  of  battery  sys¬ 
tems.  No  details  about  the  underlying  electrode  morphology 


are  involved.  While  an  obvious  weakness  of  this  approach  is 
that  it  cannot  be  used  to  design  improved  electrodes  in  an 
explicit  manner,  a  corresponding  strength  is  that  detailed  knowl¬ 
edge  of  the  cell  electrochemistry  is  not  required.  We  shall  not 
comment  further  on  the  detailed  characteristics  of  the  cell  inves¬ 
tigated  other  than  to  note  that  it  is  representative  of  batteries 
comprising  a  lithium  titanate  anode  and  cobalt  oxide  cath¬ 
ode. 

This  paper  is  organized  as  follows.  After  describing  the  exper¬ 
imental  elements  of  this  work,  we  derive  the  mathematical 
relations  employed  to  simulate  voltammetric  and  chronopo- 
tentiometric  operation  of  battery  cells.  The  following  section 
provides  a  discussion  of  results,  including  the  comparison  of 
closed-form  analytic  solutions  with  experimental  data.  We  close 
with  a  short  summary  and  a  description  of  open  questions  and 
remaining  challenges. 

2.  Experimental 

Over  the  voltage  range  of  1 .8-2.8  V,  the  prismatic  cells  main¬ 
tain  a  capacity  of  approximately  2.7  Ah.  Prior  to  each  of  the 
cyclic  voltammetry  tests,  the  cell  was  charged  at  the  C/10  rate 
(0.27  A)  to  2.8  V,  held  at  that  voltage  until  the  current  dropped 
to  less  than  50  mA,  and  subsequently  discharged  at  a  C/10  rate 
for  2,  4,  6,  or  8  h  to  reach  depths  of  discharge  (DOD)  of  20,  40, 
60,  or  80%,  respectively.  The  cell  was  placed  in  the  open-circuit 
condition  for  at  least  15  min  before  the  cyclic  voltammetry  anal¬ 
ysis. 

Scan  rates  during  cyclic  voltammetry  ranged  from  0.5  to 
5mVs-1  and  the  voltage  was  scanned  between  —15  and 
+15  mV  versus  open-circuit  voltage.  Slower  scan  rates  were  not 
employed  because  the  charge  passed  during  slower  scan  rates 
led  to  significant  changes  in  the  SOC  and  thus  the  open-circuit 
potential  Vo.  Conversely,  faster  scan  rates  were  not  employed 
as  the  5  mV  s-1  scan  rate  over  the  30  mV  potential  range  for 
the  cycle  provided  a  short  half-cycle  time  6  and  a  nearly  linear 
current-voltage  responses  (i.e.,  a  substantially  greater  than  1). 
Thus,  as  will  be  made  clear  in  the  next  section,  the  scan  rates  and 
cycle  durations  allowed  for  a  sufficiently  broad  range  in  values 
of  a  =  RC/0  to  be  examined. 

Tests  were  conducted  using  a  1287  Solartron  Electrochemi¬ 
cal  Interface  (Solartron  Analytical,  UK)  controlled  by  Corrware 
(Scribner  Associates  Inc.,  Southern  Pines,  NC).  The  cell  was 
maintained  at  25  °C;  temperature  control  was  provided  by 
a  Tenney  environment  chamber  (Model  Tenney 20,  Lunaire, 
Williamsport,  PA). 


Fig.  1.  Equivalent  circuit  for  correlating  the  cell  behavior.  Positive  currents 
denote  battery  charging. 


3.  Mathematical  relations 

This  section  contains  three  subsections.  First,  we  provide  a 
means  to  describe  the  equilibrium  potential  consistent  with  the 
behavior  of  intercalation  electrodes.  Next,  analytic  expressions 
are  derived  for  cyclic  voltammetry  of  the  cell,  and  these  expres¬ 
sions  are  used  to  simulate  the  experimental  data.  Last,  because 
the  mathematics  are  quite  similar  and  represent  an  incremental 
change  to  the  voltammetry  expressions,  we  present  solutions  for 
triangular  current  sweep  chronopotentiometry. 
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Fig.  2.  Open-circuit  potential.  The  interaction  potential  is  provided  by  Eq.  (1), 
and  Eq.  (2)  corresponds  to  the  depicted  Nernst  equation. 


3.1.  Equilibrium  potential 


Time,  x  =  f/0 

Fig.  3.  Source  term.  For  both  cyclic  current  and  potential  sources  (7  and  V, 
respectively)  the  plotted  source  term  applies,  and  the  cell  is  taken  to  be  initially 
equilibrated;  i.e.,  zero  current  for  a  prolonged  duration  and  the  cell  potential 
corresponding  to  the  equilibrium  value. 


A  more  complete  description  of  equilibrium  relations  based 
on  molecular  thermodynamics  for  intercalation  electrodes  can 
be  found  in  Refs.  [26-28].  By  considering  repulsion  between 
intercalate  species  within  the  capacity-limiting  intercalation 
electrode,  we  obtain  the  following  relation  for  the  cell  potential 
[27], 

n  Rc,T  1  —  SOC 

Vo(SOC)  =  U°--L-  In  SQC  +  Y,  ajsoa,  (1) 

7=l,oo 

where  j+  1  denotes  the  number  of  intercalates  involved  in  the 
intercalate-intercalate  interaction  and  SOC  reflects  the  frac¬ 
tional  state  of  charge.  For  our  purposes,  we  truncate  the  series 
at  7  =  6.  The  Nernst  approximation  results  from  neglecting  the 
summation  term  used  to  address  lithium-lithium  interactions: 

n  RoT  1  -  SOC 

Vo(SOC)  =  U° - In - .  (2) 

F  SOC 

The  coefficients  aj  used  to  generate  the  interaction  potential 
depicted  in  Fig.  2  are  provided  in  Table  1.  The  plot  shown  in 
Fig.  2  indicates  that  interactions  between  lithium  guests  in  the 
capacity-limiting  electrode  are  quite  weak,  as  there  is  less  than 
1 .2  kcal  mol- 1  difference  between  the  Nernst  expression  and  the 
data  at  50%  SOC. 


Table  1 

Interaction  coefficients  and  other  values  used  in  Eqs.  (1)  and  (2) 


Quantity 

Value,  V 

RT/F 

0.02568 

U° 

2.3044 

a\ 

0.8889 

a2 

-4.4742 

a3 

9.4101 

a4 

-8.7417 

a5 

3.1588 

a6 

2.3354  x  10“6 

3.2.  Cyclic  voltammetry 

The  equations  governing  the  circuit  depicted  in  Fig.  1  can  be 
written  as 

I  =  h  +  h,  v  =  Vo  +  IR  +  hRct, 

Q  &Q 

V=V0  +  IR+ /2  =  -£-  (3) 

C  at 


A  current  balance  yields  the  first  equation,  summing  the 
voltage  about  loops  including  the  interfacial  charge-transfer 
resistance  Rct  and  capacitance  C  provides  the  second  and  third 
equations  (respectively),  and  the  last  equation  relates  the  charge 
Q  on  the  capacitor  in  terms  of  the  current  passing  through  it. 
Substitution  can  be  used  to  eliminate  the  currents  /2  and  /3  and 
the  charge  Q ,  which  yields  a  single  equation  relating  the  cell 
current  I  and  potential  V: 


dl  1  /  R 

R - 1 - (  1  +  — 

d t  C  V  Rct 


_  dV  1 
dt  R&C 


(Vo  -  V). 


(4) 


It  is  helpful  to  nondimensionalize  the  problem  statement  with 
the  following  definitions: 


r  =  -, 
0 


V  = 


RC 


V- Vp 

Imax  Vo  ’ 

R 


a  = 


r  = 


Rc 


I 

I  =  — , 
vC 


Referring  to  Fig.  3,  we  can  express  the  triangular  voltage 
excitation  source  by  means  of  the  Fourier  series  [29] 


8  ~  1 
V  = - J  /  t 

77-2  Z. - J  m2 


COS 


^  n * 
n= 1,3,5,... 


nn  [  x  +  - 


Eq.  (4)  can  now  be  recast  as 


dl  ~  1  (r  ~  dV 

a— +(i+r)/=  -  -V+  — 
dr  2  V  a  dr 


The  temperature  corresponds  to  25  °C. 


(5) 
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or 


d  7  ~  r4  ^  1 

a  A  h  (1  +  r)/  — - ^  2^  “2 

dr  cy  7rz  f—-'.  nz 

77=1,3,5,... 


4^1. 
+  -  >  -  sin 

n 

n= 1,3,5,... 


127T  (  T  +  - 


COS 


mt  (  r  +  - 


(6) 


The  initial  condition  corresponds  to  an  equilibrated  system, 
V=  V0,  and 

m  =  o. 

We  have  obtained  the  following  solution  to  these  equations: 


3.3.  Triangular  current  sweep  chronopotentiometry 

Eq.  (4)  remains  as  the  governing  equation,  and  the  current 
excitation  source  is  depicted  in  Fig.  3.  The  appropriate  nondi- 
mensionalization  for  the  triangular  current  sweep  is 

t  V  -  Vo  -  I 

r  =  -,  V= - 1= - , 

0  f max  R  fmax 

RC  R 


(< a2n27t 2  +  r  +  r2){e  (O+'OAO7  cos((l /2)mt)  —  cos[ii7r(r  +  (1/2))]} 


00 

~I(t)  —  4 

n= 1,3,5,... 


— amr{e  (O+^AO7  sin((l/2)i27r)  —  sin[i27r(r  +  (1/2))]} 
an27t2[a2n27t2  +  (1  +  r)2] 


(7) 


For  long  times,  a  uniform  and  sustained  periodic  state  is 
obtained,  which  is  of  most  interest  in  terms  of  comparing  exper¬ 
iment  and  theory: 


/(periodic) 


00 

*  E  - 

71=1,3,5,... 


ocmt  sin[/t7r(r  +  (1/2))] 

—(a2n27t2  +  r  +  r2)  cos[mt(z  +  (1/2))] 


an2  tv2  [a2  n2  tv2  +  (1  +  r)z] 


(8) 


Thus,  only  the  dimensionless  voltage  and  current  differ  from 
that  of  the  cyclic  voltammetry  problem,  and  an  overbar  is 
employed  to  reflect  the  nondimensionalization  of  these  variables 
for  the  triangular  current  sweep  problem.  The  current  excitation 
source  corresponds  to 


I  =  - 


7t * 


00 

E 


n= 1,3,5,... 


1 

— ~  cos 
nz 


nn 


and  Eq.  (4)  can  now  be  written  as 


While  these  equations  are  cumbersome  because  of  the  infinite 
(albeit  convergent)  summations,  particularly  convenient  solu¬ 
tions  can  be  obtained  by  evaluating  the  summations  in  Eq.  (8)  at 
selected  times  within  the  cycle.  The  current  halfway  into  each 
cycle  (r  =  1/2)  can  be  expressed  as 

~  / 1  \  r  +  r2  +  2a  tanh((l  +  r)/ 2a) 

\2/  2a(l+r)2 

We  have  also  obtained  analytic  solutions  at  the  start  and  end 
of  each  cycle, 


7(1)  =  -7(0)  = 


sech((l  +  r)/ 2a)  —  1 
d+r)2 


(10) 


We  close  this  section  by  noting  that  relative  to  the  previously 
published  capacitor  study  [1],  which  involved  the  simple  RC 
circuit  ( Rct  infinite  in  Fig.  1,  corresponding  to  r  =  0),  we  recover 
the  same  expressions  for  Eqs.  (5)-(10)  as  0.  In  addition, 
these  analytic  expression  are  in  agreement  with  those  presented 
by  Macdonald  [30],  after  taking  into  account  the  differences  in 
the  derivation  and  notation. 


dV  dl 

a - \-  rV  =  (1  +  r)I  +  a  — 


dr 


dr 


or 


dV 

a - b  rV 

dr 


-(1+rW 


tca 


oo 

£  i 


COS 


77  =  1,3,5,... 


nn  (  r  +  - 


o 

+  a—  > 

7T 


77=1,3,5,... 


-  sm 
n 


mt  (  r  +  - 


The  initial  condition  corresponds  to  a  system  equilibrated, 
V(0)  =  0. 

We  can  derive  solutions  quite  similar  to  those  for  the  cyclic 
voltammetry  case: 


00 


V(r)  =  8  ]T 

77  =  1,3,5, 


(a2n27t2  +  r  +  r2){e  cos(mt) 

—cos[mt(z  +  (1/2))]}  +  an7t{e~^°^r  sin(ii7r)  —  sin[ii7r(r  +  (1/2))]} 


n27t2(a2n27t2  +  r2) 
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The  stationary  state  (uniform  and  sustained  periodic  state) 
corresponds  to 


V  (periodic) 


arm  sin[^7r(r  +  (1/2))] 


00 

—8  E 

n= 1,3,5,... 


+(a2n27t2  +  r  +  r2)  cos[/t7r(r  +  (1/2))] 
n27t2(a2n27t2  +  r2) 


As  in  the  case  of  cyclic  voltammetry,  the  infinite  summations 
can  be  recast  as  simple  expressions  at  various  points  within  the 
cycle, 


V(l)  =  - 


Vl} 


2a  r 

y(°)  =  72 


1  —  sech  (  — 


(s)] 


r  -\-r2  —  2a  tanh(r/2a) 


(ID 


4.  Results 


The  variation  in  the  current  response  to  the  voltage  excita¬ 
tion  is  represented  by  the  plots  shown  in  Figs.  4  and  5.  The 
two  dimensionless  groups  that  govern  the  problem,  a  =  RC/0 
and  r  =  R/Rct ,  are  varied  to  examine  their  influence.  For  r-  0, 
the  previously  published  expressions  and  plots  are  obtained  [8], 
which  have  been  helpful  in  describing  the  performance  of  super¬ 
capacitors.  As  r  is  decreased,  reflecting  increased  interfacial 
resistance  RcU  the  dimensionless  current  II vC  is  deceased,  and 
the  dimensionless  voltammogram  exhibits  increased  hysteresis. 
Conversely,  decreasing  a  leads  to  increasing  dimensionless  cur¬ 
rents  II vC  and  more  ohmic  behavior;  i.e.,  increasing  linearity  in 
the  voltammogram  and  reduced  hysteresis.  Increasing  hystere- 


Fig.  4.  Cyclic  voltammogram  for  infinite  resistance  Rct.  In  this  case,  a  series 
resistor-capacitor  circuit  prevails  (cf.  Fig.  1),  and  the  analytic  solution  presented 
in  Ref.  [8]  applies. 


Potential,  (V-V0)I(V^- V0) 


Fig.  5.  Influence  of  a  =  RC/0  and  r=R/Rct  on  the  system  response. 


sis  is  reflected  by  a  larger  magnitude  in  the  current  IlvC  at  the 
potential  (V-  Vo)l(Vmax  -  Vo)  =  0. 

The  series  of  plots  depicted  in  Figs.  6-9  provides  compar¬ 
isons  between  the  analytic  stationary- state  voltammogram  (Eq. 
(8))  and  experimental  data.  In  general,  very  good  agreement 
was  obtained,  as  the  salient  features  of  the  plotted  experimental 
voltammograms  are  reflected  by  the  accompanying  calculations. 
As  the  scan  rate  is  increased,  less  hysteresis  is  displayed,  and 
the  current-voltage  relation  becomes  linear.  Because  the  cur¬ 
rent  is  scaled  with  the  scan  rate,  the  dimensionless  current  IlvC 
declines  with  increasing  scan  rate.  As  the  SOC  increases,  the 
hysteresis  increases,  which  is  reflected  by  an  increase  in  Rct  and 
a  decrease  in  C. 
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The  trend  in  the  parameters  R ,  RcU  and  C  with  variations  in 
SOC  is  illustrated  in  Fig.  10.  The  parameters  are  effectively 
constant  for  the  constant-voltage  portion  of  the  open-circuit 
potential  curve,  and  they  change  in  a  substantially  monotonic 
manner  for  SOC  values  above  50%  SOC. 

5.  Summary  and  open  questions 

For  cyclic  voltammetry  and  triangular  current  sweep 
chronopotentiometry,  a  uniform  and  sustained  periodic  state  is 
obtained  for  batteries,  which  is  of  interest  in  terms  of  comparing 


experiment  and  theory,  and  which  admits  a  convenient  analytic 
solution  that  facilitates  parameter  extraction.  Particularly  simple 
expressions  result  for  the  response  at  the  start,  end,  and  mid¬ 
point  of  the  cycle  in  the  periodic  (stationary)  state.  This  work  is 
an  extension  of  our  earlier  analysis  of  supercapacitor  electrodes 
[8]. 

The  parameters  governing  the  cell  impedance  are  constant 
over  the  constant- voltage  portion  of  the  cell  potential,  and  vary 
monotonically  throughout  the  sloping  portion  of  the  open-circuit 
potential  (i.e.,  above  approximately  50%  state  of  charge;  cf. 
Fig.  10). 
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Fig.  7.  Experiment-theory  comparison  at  40%  state  of  charge. 


Fig.  9.  Experiment-theory  comparison  at  80%  state  of  charge. 
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In  addition  to  not  addressing  microscopic  phenomena  within 
the  cell  that  provides  a  more  comprehensive  understanding  of 
the  system,  the  primary  short  coming  of  the  method  described  in 
this  work  is  the  lack  of  low-frequency  information.  The  analytic 
expressions  are  derived  assuming  the  open-circuit  potential  Vo 
is  constant  as  the  potential  is  scanned  over  a  small  range  (30  mV 
in  this  work).  It  is  clear  from  the  inspection  of  Fig.  2  that  the 
small  variations  in  potential  can  lead  to  large  variations  in  SOC, 
particularly  for  the  flat  region  of  the  curve,  and  this  limits  the 
use  of  slow  scan  rates.  A  richer  analysis  would  include  SOC 
variations  and  the  influence  of  diffusion  processes  in  a  more 
detailed  manner,  perhaps  by  incorporating  a  Warburg  impedance 
element  as  a  starting  point  [31].  Because  the  SOC  would  change 
over  the  course  of  the  cycle,  it  is  doubtful  that  useful  analytic 
formulae  can  be  derived  for  this  case;  the  resulting  numerical 
model,  however,  should  be  easily  integrated. 

While  the  parameter  analysis  presented  in  this  work  does 
address  SOC  variations  between  experiments  (cf.  Fig.  10),  the 
circuit  parameters  will  also  change  with  temperature,  and  we 
have  only  investigated  room  temperature  variations  in  this  work. 
This  discussion  highlights  the  need  to  identify  parameter  val¬ 
ues  over  a  range  of  experimental  conditions  (particularly  over 
temperature  and  SOC  ranges);  the  very  simple  expressions 
we  derive  (Eqs.  (9)  and  (10)  for  cyclic  voltammetry  and  Eqs. 
(11)  for  triangular  current  sweep  chronopotentiometry)  should 
prove  helpful  for  extracting  parameters  over  a  broad  range  of 
temperatures. 

In  closing,  we  note  that  the  ability  of  the  simple  circuit 
depicted  in  Fig.  1  to  accurately  represent  the  experimental  cyclic- 


voltammetry  data  implies  that  the  equivalent  circuit  can  be  used 
to  represent  this  class  of  titanate-based  cells  and  that  the  equiv¬ 
alent  circuit  should  work  well  for  the  purposes  of  constructing 
model-based  state  estimators  [1-5]  for  hybrid-electric  vehicle 
applications  or  other  applications  wherein  the  state  of  the  energy 
storage  device  must  be  ascertained. 
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